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A  theoretical  estimate  is  attempted  for  the  effect  of  a  compliant  coating  on  turbulent  boundary 
layer  wall  pressure  fluctuations.  The  basic  derivation  shows  that  the  problem  reduces  to  one  of 
finding  the  distribution  in  the  wall  plaue  of  two  correlations  involving  the  wall  pressure  and  its 
normal  derivative.  Exact  expressions  arederived  for  two-dimensional  traveling  wave  pressure/ve¬ 
locity  admittances  of  an  isotropic  elastic  coating.  These  admittances  are  combined  with  some 
reasonable  assumptions  about  the  form  of  the  pressure  cross  spectral  density  to  yield  approximate 
expressions  for  the  two  desired  pressure/derivative  correlations.  Finally,  two  surface  integrals  of 
these  correlations  result  in  the  wall  pressure  function  in  the  presence  of  the  compliant  boundary. 
The  calculations  indicate  that  the  compliant  wall  increases  the  mean  square  wall  pressure  at  low 
speeds  and  decreases  the  pressure  fluctuations  at  high  speeds.  Unfortunately,  the  reduction  at 
high  speeds  probably  cannot  be  achieved  in  practice  because  of  the  related  mechanical  problem 
of  static  divergence  of  the  coatiLg . 
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A  THEORETICAL  ESTIMATE  OF  TURBULENT  WALL  PRESSURE 
FLUCTUATIONS  ON  A  COMPLIANT  BOUNDARY 


INTRODUCTION 

$ 

The  problem  of  boundary  layer  behavior  in  the  presence  of  a  compliant 
boundary  was  brought  to  attention  by  the  intriguing  experiments  of  M.  O. 

Kramer, 1 » 2>  *  who  demonstrated  reduced  drag  on  underwater  towed  bodies 
covered  by  flexible  coatings.  Kramer  intuitively  ascribed  the  drag  reduction  to 
a  transition  delay  provided  by  the  dissipation  of  the  compliant  boundary.  Sub¬ 
sequent  theoretical  studies,  using  an  extension  of  the  Orr-Sommerfeld  stability 
equation,  verified  the  transition  delay  but  indicated  that  the  cause  was  not  added 
dissipation  but  rather  a  profound  modification  of  the  disturbance -wave  structure 
of  the  flow.  These  linearized  laminar  stability  studies  were  begun  by  Boggs 
and  Tokita4  and  later  improved  and  extended  by  Benjamin,  B»  *  Nonv  xler, 7 
Linebarger,8  Landahl , 9  Hains, 10  and  Kaplan. 11  Kaplan’s  thesis  summarizes 
the  previous  work  and  contains  extensive  numerical  stability  calculations  for  a 
variety  of  model  compliant  boundaries. 

Although  the  linearized  theories  definitely  predict  transition  delay,  the  ex¬ 
periments  which  have  followed  Kramer's  pioneering  work  have  met  only  mixed 
success.  The  measurements  of  VonWinkle, 12  Boggs  and  Frey, 13  and  Laufer 
and  Maestrello14  do  not  yield  a  definitive  interpretation,  largely  because  the 
transition  delay,  if  it  exists,  is  difficult  to  separate  from  whatever  effect  the 
coating  might  have  on  the  fully  turbulent  region.  It  is  hoped  that  measurements 
in  the  new  water  tunnel  at  the  Underwater  Sound  Laboratory  might  clarify  the 
matter. 

The  fully  turbulent  boundary  layer  in  the  presence  of  a  compliant  surface 
presents  a  formidable  theoretical  challenge .  The  sound  radiated  by  turbulence 
near  a  flexible  boundary  has  been  estimated  by  Ffowcs  Williams  and  Lyon, 15 
and  the  Reynolds  stress  very  near  the  surface  has  been  studied  by  Ffowcs 
Williams. 16  To  the  authors'  knowledge,  there  has  been  no  theoretical  study  of 
turbulent  wall  pressure  fluctuations  at  a  compliant  boundary,  and  it  is  the  intent 
of  the  present  report  to  present  such  a  theory.  The  general  development  here¬ 
in  can  easily  be  sifted  for  qualitative  information.  However,  to  produce  quan¬ 
titative  results,  the  analysis  has  resorted  to  a  series  of  approximations  which 
probably  reduce  the  final  calculations  to  the  status  of  a  fairly  crude  estimate. 
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BASIC  ANALYSIS 

The  basic  problem  of  turbulent,  boundary  layer,  wall  pressure  fluctuations 
has  been  the  subject  of  intense  experimental  study  over  the  past  decade  for  the 
particular  case  of  a  rigid  wall.  Corcos 17  gives  an  excellent  review  of  the  many 
published  measurements  of  fluctuating  pressure  on  a  rigid  wall.  Probably  the 
most  accurate  of  these  measurements  are  those  of  Bakewell  et  al. 18  and  Will- 
marth  and  Wooldridge, 19  and  it  will  be  necessary  to  use  these  data  in  the  anal¬ 
ysis  which  follows.  To  the  authors'  knowledge,  no  measurements  of  fluctuating 
pressure  at  a  compliant  wall  have  appeared  in  the  open  literature. 

Although  measurements  abound,  theoretical  work  on  turbulent  wall  pres¬ 
sure  is  lacldng.  Based  on  an  approach  suggested  by  Gardner,  20  > 21  a  complete, 
though  rather  approximate,  theory  of  the  space -time  distribution  of  rigid  wall 
pressure  has  been  given  by  White. 22  White's  analysis  indicates  that  the  sta¬ 
tistical  properties  of  pressure  at  a  rigid  wall  are  primarily  affected  by  the 
shape  of  the  mean  velocity  profile  in  the  boundary  layer.  In  particular,  White's 
results  predict  that  the  longitudinal  space  correlation  is  affected  significantly 
by  the  mean  velocity  profile,  while  the  lateral  correlation,  the  power  spectrum, 
and  the  convection  speeds  are  affected  very  little.  Actual  measurements  seem 
insensitive  to  profile  shape,  a  phenomenon  Corcos17  calls  "space-time  simi¬ 
larity,"  although  recent  data  by  Schloemer 22  for  pressure  gradients  indicate 
some  profile  effect,  particularly  on  the  power  spectrum.  This  apparent  overall 
insensitivity  of  rigid  wall  data  is  exploited  in  the  present  analysis. 

The  fluctuating  pressure  p  may  be  calculated  in  principle  for  incompres¬ 
sible  flow  of  a  Newtonian  fluid  by  taking  the  divergence  of  the  Navier-Stokes 
equations,  yielding  the  Poisson  equation 


V'-/?  =  — p  S(X,/) 


(1) 


where  the  function  S  is  a  complicated  combination  of  velocity  derivatives  and 
t  is  time.  The  actual  form  of  S  is  given  by  White22  and  is  not  important  in 
the  present  study.  The  position  vector  X  has  coordinates  (x1,x2,x.1)  which 
are  sketched  in  Fig.  1.  The  freestream  flows  in  the  x,  direction. 

The  formal  solution  of  Eq.  (1)  for  pressure  at  the  wall  (x2  =  0)  is 
given  by  Green's  function  integral  solution: 
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Note  that  the  second  integral  in  Eq.  (2)  requires  knowledge  of  a  boundary  con¬ 
dition  in  the  form  of  the  normal  derivative  of  p  at  the  wall.  For  a  rigid  wall, 
this  derivative  is  negligibly  small ,  by  analogy  with  boundary  layer  theory,  as 
Kraichnan24  has  shown.  Thus,  the  rigid  wall  pressure  is  given  simply  by  the 
first  integral,  which  involves  only  the  source  term  S.  If  we  accept  the  exper¬ 
imental  evidence  that  the  wall -pressure  correlation  is  insensitive  to  the  form  of 
S,  then  the  effect  of  a  compliant  wall  must  be  primarily  due  to  changes  in  the 
boundary  condition  on  p.  It  is  the  purpose  of  this  report  to  investigate  how  the 
compliant  surface  might  affect  the  normal  derivative  of  p  at  the  wall,  so  that 
the  mean  value  of  the  second  integral  in  Eq.  (2)  might  be  evaluated,  at  least 
approximately. 

To  shorten  the  expressions  which  follow,  let  P  denote  wall  pressure  p7(. 
and  let  us  rewrite  Eq.  (2)  with  the  following  tighter  notation  for  the  double 
and  triple  integrals: 


j_  r  /  dP  \  dz  P  r  sdz 
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If  we  define  the  wall -pres  sure,  space -time  correlation  by  the  relation 


R(X.X',t,i')  =  P(X,t)P(X',i')  ,  (4) 


where  the  overbar  denotes  the  time  average  in  the  statistically  stationary  sense, 
then,  by  substituting  into  Eq.  (3)  and  performing  the  time  average  underneath 
the  integral  signs,  we  obtain 


3 


where  the  integrals  involving  P  are  to  be  evaluated  in  the  planes  (z2fz2'  =  0) 
and  the  integrals  involving  S  are  evaluated  in  the  infinite  half  space  above 
this  plane  of  the  wall.  After  inspection,  we  find  that  the  first  two  integrals  on 
the  left-hand  side  are  identical  because  of  symmetry  in  a  plane. 

It  should  be  noted  that  the  right-hand  side  involves  the  source  terms  S 
which  occur  in  the  boundary  layer  flow  past  whatever  type  boundary  is  under 
study.  That  is,  if  we  seek  to  use  Eq.  (5)  to  calculate  R  for  a  compliant 
surface,  then  S  should  be  the  source  function  for  flow  past  a  compliant  sur¬ 
face.  It  is  at  this  point  that  we  use  the  experimental  insensitivity  of  the  source 
function,  previously  discussed,  to  postulate  that  the  right-hand  side  of  Eq.  (5) 
is  essentially  identical  to  R0,  the  pressure  correlation  in  the  presence  of  a 
rigid  wall.  This  assumption,  although  reasonable,  cannot  be  verified  until  data 
are  available  for  mean  and  fluctuating  velocities  in  the  boundary  layer  past  a 
compliant  surface.  Apparently  Professor  J.  Lumley  at  the  Pennsylvania  State 
University  is  presently  making  such  measurements.  The  question  is  also  being 
examined  theoretically  at  present  by  the  second  author  as  a  thesis  for  the  Uni¬ 
versity  of  Rhode  Island. 

Combining  the  first  two  integrals  in  Eq.  (5)  and  utilizing  the  assumption 
that  the  source  integral  is  equal  to  R0,  one  obtains  the  following  basic  rela¬ 
tion  for  calculating  the  wall  pressure  at  a  compliant  surface: 


dP  dP' 
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(6) 


At  first  glance,  Eq.  (6)  might  appear  to  predict  that  the  compliant  wall  cor¬ 
relation  R  is  always  less  than  the  rigid  wall  value  R0.  However,  we  shall 
see  that  the  first  integral  is  usually  negative,  while  the  second  integral  is  posi¬ 
tive,  with  the  result  that  the  effect  on  R  is  rather  mixed. 
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TRAVELING  WAVE  ADMITTANCES 


Since  R0  is  known  from  experiment,  the  evaluation  of  R  from  Eq.  (6) 
can  be  accomplished  if  the  correlations  involving  P  and  its  normal  derivative 
can  be  estimated.  To  do  this,  we  must  investigate  the  properties  of  an  idealized 
elastic  coating.  All  available  experiments  indicate  that  turbulent  boundary 
layer  fluctuating  pressures  have  approximately  the  form  of  traveling  waves 
moving  in  the  X;  direction  with  a  convection  velocity  Ue  somewhat  less 
than  the  ireestream  velocity  Ua  .  Naturally,  there  is  a  certain  amount  of  con¬ 
vective  incoherence,  since  the  pressure  waves  as  they  move  downstream  are 
undergoing  continuous  decay  and  regeneration.  No  attempt  will  be  made  here 
to  reproduce  this  effect;  that  is,  the  fluctuating  pressures  will  be  treated  as  a 
simple  summation  of  traveling  waves  of  different  frequencies.  A  second  diffi¬ 
culty  is  that  the  actual  turbulent  pressures  are  not  plane  waves  but  instead  have 
some  unknown  variable  shape  in  the  lateral  (x3)  direction.  This  analysis 
will  treat  the  case  of  plane  waves  and  then  attempt  belatedly  to  introduce  a 
three-dimensional  effect  by  use  of  the  measured  lateral  spectra  of  wall  pres¬ 
sure. 

Consider  a  compliant  coating  of  thickness  h,  backed  up  by  a  rigid  under¬ 
surface,  as  shown  in  Fig.  1.  Several  studies  have  been  made  of  the  response 
of  such  a  coating  to  a  plane  traveling  wave  for  a  Hookean  isotropic  coating,  as¬ 
suming  small  strains.  The  analytical  results  are  in  the  form  of  traveling  wave 
admittances,  which  are  amplitude  ratios  of  coating  velocity  to  traveling  wave 
pressure.  Following  a  suggestion  of  Nonweiler, 7  Kaplan"  calculated  admit¬ 
tances  by  assuming  a  condition  of  plane  stress  in  the  coating,  while  Tokita  and 
Boggs25  gave  admittances  for  the  case  of  plane  strain.  Both  Kaplan's  results 
and  those  of  Tokita  and  Boggs  contain  algebraic  errors  which,  hopefully,  have 
been  eliminated  in  the  present  report.  Also,  Tokita  and  Boggs,  by  expanding  in 
a  series  and  truncating,  gave  approximate  admittances  (equation  7b9  of  Ref. 

25)  which  they  later  used  in  a  study  of  coating  stability. 26  However,  numerical 
evaluation  of  exact  admittance  formulas  shows  that  these  approximations  are 
valid  only  for  a  small  range  of  frequencies  and  hence  will  not  be  used  here. 

As  is  usual  in  elasticity  theories,  there  is  no  great  difference  between  the 
plane  stress  admittances  and  the  plane  strain  results.  Let  us  reproduce  a  plane 
stress  analysis,  similar  to  that  of  KapTan,  "  comparing  the  final  results  ob¬ 
tained  with  those  of  Tokita  and  Boggs25  for  plane  strain. 

From  Fig.  1,  a  streamwise  plane  traveling  pressure  wave  would  have  the 
following  complex  form  for  any  given  frequency  u>  : 


in  (jc,  —  Urt) 


(7) 


P  —  Pi,  e 


where  Ur  is  the  convection  speed  and  «  =  io/Uf  is  the  wave  number.  Let 
t  be  the  coating  displacement  in  the  x,  direction  and  >|  be  the  displace¬ 
ment  in  the  x2  direction.  Let  the  elastic  coating  have  shear  modulus  G, 
Poisson’s  ratio  n,  and  Young's  modulus  E  =  2G(1  +  n).  Then  the  equations 
of  elasticity  for  plane  stress  are 
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Equation  (8)  contains  eight  linear  algebraic  and  differential  equations  in  the 
eight  variables  on,  oaa,  o1:!,  e, ,,  e*,,  e , »|,  and  c.  The  stresses  and  strains 
may  easily  be  eliminated  in  favor  of  the  two  displacements  for  which  boundary 
conditions  are  known  at  the  upper  and  lower  surface  of  the  coating.  Since  the 
system  is  linear  and  the  driving  force  is  a  traveling  wave  (Eq.  (7) ) ,  it  follows 
that  the  resulting  displacements  must  also  be  traveling  waves  with  amplitudes 
which  vary  through  the  thickness.  Hence  we  postulate  that 
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The  exponential  expressions  will  cancel  properly  from  the  equations  of  motion, 
leaving  a  single  ordinary  fourth  order  linear  differential  equation  in  i|„ 
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where  the  constants  rt  and  r,  are  related  to  the  ratio  of  the  convection  speed 
Ur  to  the  coating  shear  wave  speed  C„  =  \  g as  follows: 


r,-=  !  -  '/2  (I  -f0  U,-/ c;- 
r,.-  =  I  -  u,  -/c:- . 


(11) 
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The  primes  in  Eq.  (10)  indicate  differentiation  with  respect  to  x2.  An  equa¬ 
tion  identical  to  Eq.  (10)  holds  for  the  other  displacement,  So. 

The  general  solution  of  Eq.  (10)  is 

i|„  =  A,  sinh  (dr,.v_.)  -(-/T.cosh  (tu\x..)  -f  d:isinh  ( ar.x-< )  -f-/j4cosh  (a r-.x.)  (12) 


where  A, — A 4  are  constants.  Assuming  that  the  coating  is  securely  bonded 
+o  the  rigid  under  structure,  the  boundary  conditions  at  the  lower  surface  state 
that  the  displacements  must  vanish: 


lA-h)  =  I),,  (-h)  =0. 


(13) 


At  the  upper  s  urface,  the  vertical  normal  stress  in  the  coating  must  equal  the 
traveling  wave  pressure: 
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A  simple  and  realistic  fourth  boundary  condition  is  achieved  by  setting  the  shear 
stress  equal  to  zero  at  the  upper  surface: 


+  _9l  =  o  at  x,  =  0  . 
3*"  3*i 


(15) 


Actually,  the  shear  stress  at  the  upper  surface  does  not  vanish  but  instead  must 
equal  the  fluid  shear  stress  in  the  boundary  layer  a't  the  wall,  rw.  However, 
tn-  is  small,  and  the  traction  of  t,f  assigned  to  any  given  traveling  wave 


must  be  very  small  indeed.  Hence  this  refinement  is  not  considered  to  be  nec¬ 
essary  for  an  accurate  calculation  of  the  admittance  of  the  coating. 

Equations  (13),  (14),  and  (15)  are  sufficient  to  define  unique  values 

of  the  four  constants  A,  in  Eq.  (12).  Since  we  will  ultimately  be  concerned 
with  velocities  at  the  upper  surface,  it  is  convenient  to  give  the  solutions  in 
terms  of  the  admittances  Y„  and  Y, ,  defined  as  follows: 


(16) 


The  negative  sign  in  the  definition  of  Y„  is  traditional.  Unfortunately,  the 
admittances,  although  defined  as  ratios,  are  not  dimensionless.  It  will  be  con¬ 
venient  to  use  the  dimensionless  group  ( P ,  C.  Y)  ,  which  is  a  function  of  the 

dimensionless  parameters  r  x ,  r  2 ,  and  ( a  h) .  The  normal  admittance  is 
given  by  the  expression 


/' 


C  Y  — 

a  1  ii  — 


-/>,  ((7,/C,):i 
2(1  +  r/)  A ,  4-  4  r,  r-,  A;t 


(17) 


The  constants  A2  and  A4  have  been  incorporated  into  Eq.  (17),  but  the 
expressions  for  A,  and  A:j  are  rather  lengthy.  If  we  adopt  the  short  no¬ 
tation 

C j  =  cosh  (ar,/i) 

for  /  =  1,2,  (18) 

Ti=  tanh  (ttr,h) 


then  A,  and  A:i  may  be  written  as  follows: 


A  n 


-  f i  r.  7,  T.  -  '/2  ( 1  +  r.r)  ( 1  -77)  C,/C2 
r  .  -  r,  r  .  7', 

(I  —  T.)  C./C,  -'/2(l-f  r-y) 


A  i  — 


T , 


(19) 
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In  the  limit  as  (ah)  approaches  infinity,  T,  and  T2  approach  unity. 

By  inspection,  we  see  that  A3  will  approach  unity  and  A,  becomes 
~  Vi  (1+r*).  From  Eq.  (17),  the  admittance  will  approach  the  limiting  value 


C.  Y„  ( oc  ) 


-;  r,  ( UJC 
4  rt  r  -j  — -  ( 1  +  r )  - 


(20) 


The  denominator  of  Eq.  (20)  becomes  zero,  giving  infinite  Y„ ,  at  a  speed 
ratio  (Uf/C.)  varying  from  0.874  for  n  =0  to  a  value  of  0.933  for  n  =  0.5. 
For  speeds  less  than  this  critical  value,  the  denominator  is  positive,  and  Eq. 
(17)  predicts  in  general  that  Yn  will  be  a  pure  negative  imaginary  quantity 
for  any  subcritical  speed. 

Equation  (17)  applies  for  a  coating  which  is  perfectly  bonded  to  the  un¬ 
dersurface  ,  i.e.,  it  satisfies  Eq.  (13).  If  one  relaxes  this  condition,  a  much 
simpler  expression  for  Y„  results,  as  shown  by  Kaplan. 11  Instead  of  being 
bonded,  we  could  postulate  that  the  coating  slides  without  shear  along  the  lower 
surface,  satisfying  the  following  conditions: 

=0.  (21) 

dx-j 


The  use  of  Eq.  (21)  instead  of  Eq.  (13)  gives  a  much  simpler  normal  ad¬ 
mittance,  which  we  term  the  "shearless"  coating  result: 


ft.  C.  Y„  (shearless)  = 


—  i  r i  (U,./C _ 

4  r,  r,/T-,  -  (1  +  r2-’)-yr,  ’ 


(22) 


Clearly,  the  limit  of  Eq.  (22)  as  (ah)  approaches  infinity  is  identical  to 
Eq.  (20)  for  the  bonded  coating.  In  general,  for  a  given  subcritical  speed, 
there  is  no  great  difference  between  the  bonded  and  the  shearless  coating  over 
the  entire  frequency  range,  as  Fig.  2  shows,  using  m  =  0.5  as  an  approxi¬ 
mate  value  for  natural  rubber.  As  Fig.  2  indicates,  the  shearless  admittances 
at  low  frequencies  are  about  twenty  per  cent  higher  than  the  bonded  values. 

The  high  frequency  asym  !otes  are  identical. 

The  tangential  admittance  Yt  as  defined  inEq.  (16)  may  also  be  calculated. 
The  result  for  the  bonded  coating  is 


9 


(23) 


-  (  U  r/Cn  )  An  A  ]) 

2  /*]  r-j  A  3  -f"  A  \  ( ]  -j"  r->- ) 


where  A,  and  A:i  are  again  defined  by  Eq.  (19).  Once  again  a  simpler 
expression  results  for  the  "shearless”  coating: 


f>.C»  Y i  (sheariess) 


-  (Ur/C,)  (2  r,  r,  —  ( 1  +  r,-)  TJT, ) 
4  r,  rs  —  ( 1  +  rJy-  T-./T, 


(24) 


As  before,  the  admittance  is  slightly  larger  for  the  shearless  coating  as -com¬ 
pared  to  the  bonded  value.  Figure  3  compares  the  tangential  and  normal  admit¬ 
tances  for  ;he  shearless  coating  for  n  =  0.5.  The  asymptotic  values  of  Y, 
are  roughly  one -half  of  the  asymptotic  magnitude  of  Y„  for  the  same  speed 
ratio.  The  tangential  admittance  suffers  a  singularity  at  the  same  "critical" 
speed  ratio  as  Y„,  as  listed  in  Table  1. 


Table  1 

MINIMUM  SPEED  RATIO  FOR  AN  ADMITTANCE  SINGULARITY  TO  OCCUR 


Poisson's  Ratio 

Minimum  (Ur/C,) 

0.0 

0.8740 

0.1 

0.8913 

0.2 

0.9052 

0.3 

0.9162 

0.4 

0.9252 

0.5 

0.9325 
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As  we  shall  see  in  the  next  section,  the  evaluation  of  Eq.  (6)  for  a  rea¬ 
sonably  thick  coating  (of  the  order  of  the  boundary  layer  thickness)  depends 
only  upon  the  asymptotic  values  of  Yn.  Figure  4  shows  the  magnitude  of  these 
asymptotic  admittances  for  subcritical  speeds.  Note  that,  for  low  speeds,  the 
asymptotic  admittances  vary  linearly  with  speed  ratio. 

ADMITTANCE  SOLUTION  FOR  PLANE  STRAIN 

The  previous  theoretical  admittances,  Eqs.  (17)  through  (24),  are  de¬ 
rived  for  the  assumption  of  plane  stress  (zero  stress  in  the  z  direction). 

The  analogous  solution  for  plane  strain  (zero  z  displacement)  was  given  by 
Tokita  and  Boggs,25  following  a  somewhat  more  complicated  analysis,  using 
the  three-dimensional  wave  equation  which  results  from  the  definition  of  the  so- 
called  "displacement  potentials."  The  boundary  conditions  used  were  Eqs.  (13) 
and  (15),  that  is,  a  tightly  bonded  coating.  An  exact  expression  for  the  admit¬ 
tance  was  not  given  but  can  easily  be  calculated  from  equation  (7b6)  of  Ref.  25. 
The  parameter  r2  is  the  same  as  for  plane  stress,  but  the  quantity  r,  is 
slightly  differ  ent .  That  i  s , 


rr  =  1  -  u;-/c;- 

r,*-'  =  1  -  '/2  U,-/C'- , 

(1  —  n) 


(25) 


where  the  asterisk  is  included  in  r*  as  a  reminder  that  it  is  the  plane  strain 
value.  Using  this  notation,  the  exact  expression  for  the  normal  admittance  for 
plane  strain  in  a  bonded  coating  is 


P‘ i 


C,  Yn  = 


in*  ( Ur/C.Y '  (7V 


r-.  r_.) 


(26) 


4 rt*  r-,  (1  —  r,*  r.  T,*  T-,)  —  (1  +  r.-)-  (r,*  r,  —  T\*  T.)  —  4 r,*  r-,  ( 1  +  r/)/C t*  Cj 


where  T,  and  C;  are  as  defined  in  Eqs.  (18).  Figure  5  compares  values 
of  Y„  from  Eq.  (26)  to  equivalent  values  for  the  plane  stress  case,  Eq. 

(17).  For  a  given  speed  ratio,  the  plane  strain  admittance  is  somewhat  smaller 
and  has  a  lower  asymptote.  The  high  frequency  asymptote  of  Eq.  (26)  is 


C,  Y„  (  x  ) 


-in*  (U,/C,)* 

4  r,*  r-,  —  ( \  +  r-f)- 


(27) 


which  is  identical  in  form  to  Eq.  (20)  for  the  plane  stress  case.  Table  2  gives 
a  comparison  of  these  asymptotic  values. 
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Table  2 

ASYMPTOTIC  NORMAL  ADMITTANCES  (,.C.  Yn  (»)//)  FOR  n  =  0.5 


U./C. 

PLANE  STRESS 

PLANE  STRAIN 

0.2 

-  0.1374 

-0.1031 

0.4 

-  0.3038 

-0.2282 

0.6 

-  0.5676 

-0.4232 

0.7 

-  0.8118 

-0.5950 

0.8 

-  1.3406 

-0.9302 

0.85 

-  2.0189 

-1.2926 

0.9 

-  4.6190 

-2.2262 

0.91 

-  6.4855 

-2.6433 

0.92 

-11.3029 

-3.2873 

0.93 

-54.1007 

-4.4247 

Although  the  plane  strain  values  in  Table  2  are  substantially  smaller  in  mag¬ 
nitude,  we  shall  see  that  this  has  no  great  effect  on  the  wall  pressure  analysis 
which  follows.  However,  since  a  practical  coating  construction  would  probably 
be  constrained  in  a  manner  somewhere  in  between  these  two  extremes,  one  can 
look  upon  Table  2  as  a  measure  of  the  uncertainty  involved  in  a  theoretical  esti¬ 
mate  of  the  actual  coating  response  to  traveling  waves. 

As  mentioned  before,  Ref.  25  did  not  attempt  to  calculate  the  exact  plane 
strain  admittance  as  given  by  Eq.  (26).  Instead,  Tokita  and  Boggs  appro  i- 
mated  the  hyperbolic  functions  by  the  first  two  terms  of  their  Taylor  series  ex¬ 
pansions.  The  result  was  an  approximate  admittance  expression  (equation  (7b9) 
of  their  report).  In  the  present  notation,  this  approximation  is  written  as 
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where 


_ -  2  a:<  h  ' _ 

*  C '  Y”  “  3  (a-  h-  -  2)  ( U.-/C;-  -k-)' 


(28) 


_ 4  a-  h- _ 

a-  h-  (3  —  4|i)  —4(1  —  |x) 


Figure  6  compares  the  exact  admittance  from  Eq.  (26)  with  the  approximate 
value,  Eq.  (28),  for  the  case  ( \5c/C ,)  =  0.5  and  n  =0.5.  It  is  seen  that 
Eq.  (28)  is  accurate  only  for  a  small  intermediate  frequency  range.  Note  that 
Eq.  (28)  fails  to  predict  a  constant  asymptotic  admittance  at  high  frequencies. 
The  approximate  admittance,  although  apparently  rather  crude,  was  used  by 
Tokita  and  Boggs  in  Ref.  26  to  predict  the  mechanical  stability  (static  diver¬ 
gence  and  flutter)  of  a  compliant  coating.  Since  their  calculations  were  rather 
complex  and  also  involved  further  approximations,  it  is  not  clear  exactly  what 
quantitative  effect  the  error  inherent  in  Eq.  (28)  would  introduce  into  the  re¬ 
sults  of  Ref.  26. 

Finally,  we  may  note  that,  for  "supercritical"  speeds  (greater  than  those 
in  Table  1),  all  of  the  admittance  expressions  possess  multiple  singularities. 
Since  Ref.  26  predicts  a  statically  unstable  coating  at  such  speeds,  no  super¬ 
critical  calculations  were  made  in  this  report. 

PRESSURE  DERIVATIVE  CORRELATIONS  AT  THE  COMPLIANT  WALL 

The  chief  result  of  the  basic  analysis  section  of  this  report  was  to  show 
that  the  problem  of  estimating  compliant  surface  pressure  fluctuations  reduces 
approximately  to  the  evaluation  of  Eq.  (6).  The  first  integral  in  Eq.  (6)  can¬ 
not  be  evaluated  until  we  know  the  distribution  of  the  correlation  function 
P  O P')/dx-/  in  the  plane  of  the  wall.  The  second  integral  requires  knowledge  of 
the  correlation ( gp  gp ) /( g*,  gjr^in  the  wall  plane.  It  is  the  purpose  of  this  sec¬ 
tion  to  show  that  these  correlations  can  be  reasonably  approximated,  using  the 
traveling  wave  admittance  approach. 

It  is  obviously  necessary  to  the  admittance  approach  that  we  assume  that  the 
turbulent  pressure  disturbances  are  in  the  form  of  a  superposition  of  many 
small  traveling  waves  having  different  amplitudes  and  frequencies.  This  is 
certainly  not  true  on  an  instantaneous  basis.  That  is,  turbulent  pressure  fluc¬ 
tuations  suffer  by  nature  a  convective  incoherence.  The  disturbances  are  con¬ 
stantly  decaying  and  being  regenerated  as  they  move  downstream  with  a  con¬ 
stantly  changing  convection  speed.  It  is  only  on  a  time-averaged  basis  that  the 
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pressure  simulates  in  any  way  a  sum  of  traveling  waves.  Some  evidence  of 
convective  incoherence  persists  even  on  a  time -average  basis.  For  example, 
the  convection  speed  Ur  is  not  truly  constant  but  instead  varies  with  the  fre¬ 
quency  and  with  the  spacing  between  correlated  points.  Also,  the  sharp  de¬ 
crease  in  the  pressure  correlation  with  lateral  spacing  indicates  that  the  as¬ 
sumption  of  plane  traveling  waves  is  not  very  accurate,  even  on  the  average. 
However,  it  is  fortuitous  that  these  deviations  from  ideal  traveling  wave 
behavior  do  not  have  a  strong  effect  on  the  behavior  of  a  compliant  boundary, 
because,  as  the  calculations  will  show,  the  compliant  wall  responds  in  an  ex¬ 
tremely  localized  fashion  to  the  pressure  disturbances.  That  is,  the  correla¬ 
tions  needed  in  Eq.  (6)  drop  off  so  rapidly  with  distance  that  their  effect  on  the 
calculation  of  R  in  Eq.  (6)  is  confined  to  a  small  local  region  whose  diam¬ 
eter  is  less  than  a  boundary  layer  thickness.  Under  these  conditions,  the  con¬ 
vective  incoherence,  which  occurs  on  a  somewhat  larger  scale,  does  not  cause 
any  great  error  in  the  analysis. 

The  normal  derivative  of  P  is  related  to  the  velocity  components  through 
the  normal  component  of  the  Navier-Stokes  equations: 


1 2L 

t>  d-X-J 


du-j 

d 


M  V-  <«a)  • 


(29) 


For  a  xigid  wall,  u,  and  u2  both  vanish  at  the  wall,  leaving  only  the  viscous 
term  on  the  right-hand  side  of  Eq.  (29).  As  mentioned  before,  Kraichnan 24 
showed  this  viscous  term  to  be  negligibly  small  for  a  rigid  wall.  However,  for 
a  compliant  wall,  none  of  the  velocity  terms  in  Eq.  (29)  vanish,  and  care  must 
be  taken  to  ascertain  their  magnitude.  The  no-slip  condition  should  still  be 
valid,  so  that  the  fluid  velocities  at  the  wall  must  equal  the  surface  velocities  in 
the  coating,  which  in  turn  are  related  through  the  admittance  functions  to  the 
fluctuating  wall  pressure.  The  use  of  coating  velocity  instead  of  fluid  velocity 
allows  us  to  ignore  the  interplay  between  the  fluid's  mean  and  fluctuating  veloci¬ 
ties  -  an  interplay  which  has  caused  erroneous  results  in  stability  studies,  e.g., 
Ref.  4. 

To  evaluate  the  terms  in  Eq.  (29),  consider  first  a  single  traveling  wave 
of  amplitude  P„ .  Using  the  admittance  concept,  one  can  calculate  the  ampli¬ 
tude  of  the  normal  acceleration  at  the  coating  surface : 


djh 

Y„ 

d' 

B  n 

(30) 


14 


•  I 


In  a  similar  manner  we  obtain  an  estimate  of  the  first  convective  term  in  Eq. 
(29): 


Y„  Y, 

9*1 

(31) 


Let  us  denote  the  dimensionless  tangential  admittance  by  the  symbol 
Yt*  =  p ,  C,  Yf.  Then  the  ratio  of  these  two  terms  is,  approximately: 


9 Ms 
9*i 


9«a 


=  Yt* 


P„ 


PH 


C,  Ur 


(32) 


From  Fig.  4,  Yt  is  less  than  5  for  Uf/C,  less  than  0.8,  while  the  dimen¬ 
sionless  pressure  amplitude  P „/(/>„  C,  Ur)  is  much  smaller  than  unity.  Then, 
for  a  single  traveling  wave,  the  first  convective  term  is  negligibly  small  com¬ 
pared  to  the  local  normal  acceleration.  Then,  by  superimposing  a  large  number 
of  traveling  waves,  one  arrives  at  the  root -mean-square  approximation: 


All  available  measurements  indicate  that  the  root -mean-square  turbulent  pres¬ 
sure  Prm,  =  0.003  p  U*  >  where  p  is  the  fluid  density  and  Uffi  is 

the  freestream  velocity.  Thus,  the  dimensionless  pressure  in  Eq.  (33)  is  a 
very  small  fraction  for  subcritical  speeds,  making  the  first  convective  term 
negligible.  A  similar  comparison  of  the  second  convective  term  to  the  local 
acceleration  yields  exactly  the  same  order  of  magnitude  estimate  as  that  of  Eq. 
(33),  so  that  this  term  is  also  quite  small.  Finally,  the  ratio  of  the  viscous 
term  in  Eq.  (29)  to  the  local  acceleration  is  found  to  be  of  order  c/,  the 
local  skin  friction  coefficient.  Since  cf  for  a  turbulent  boundary  layer  is  ap¬ 
proximately  0.005  or  less,  the  viscous  term  is  also  negligible.  Clearly,  then, 
the  pressure  normal  derivative  in  Eq.  (29)  is  dominated  by  the  local  normal 
acceleration,  and  an  accurate  estimate  to  the  first  of  the  two  desired  pressure 
correlations  is: 
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p 


(34) 


'dPj  ^  _  pdu£ 
P  dt' 


To  evaluate  Eq. 
tion  between  P 


(34),  we  note  that,  for  a  single  traveling  wave,  the  correla- 
and  the  normal  acceleration  would  be 


Po 


dU‘2 

3t 


=  /  Y„  o>  p,r 


(35) 


To  generalize  this  expression  to  a  :omplete  distribution  of  traveling  waves,  we 
make  use  of  the  space -frequency  correlation  r  of  the  wall  pressure,  which 
is  the  Fourier  transform  of  the  space -time  correlation  R  defined  in  Eq.  (4); 

R  (1  *),  to)  =  fX  T  (£,  t|,  eo)  e' dio  ,  (36) 

where  1  ,  n  »  and  t0  are  the  longitudinal  separation,  the  lateral  sepa¬ 
ration,  and  the  time  delay  between  the  two  correlated  boundary  points,  respec¬ 
tively.  Utilizing  this  function  V  to  generalize  Eq.  (35),  we  obtain  the  fol¬ 
lowing  expression  for  the  first  desired  pressure  correlation  function: 


P 


dP' 

3*2 


re" . da), 


(37) 


where  Yn  is  taken  to  be  frequency  dependent  as  given  by  Eq.  (17),  for  ex¬ 
ample,  for  a  bonded  coating.  Equation  (37)  is  not  an  exact  representation  un¬ 
less  Yn  is  given  spatial  properties  to  account  for  the  fact  that  turbulent 
pressure  disturbances  are  not  purely  plane  waves.  The  authors  have  not  at¬ 
tempted  to  introduce  such  a  sophisticated  admittance  function  into  this  analysis, 
arguing  in  the  previous  paragraphs  that  the  ’’localized"  behavior  of  the  coating 
makes  a  spatially  distributed  admittance  unnecessary. 

Extensive  measurements  are  available  for  the  frequency  correlation  F0 
for  the  case  of  a  rigid  wall.  For  zero  separation,  I’0  reduces  to  the  power 
spectrum  <*>«  =  r0(o,o,  to).  The  data  of  Ref.  18  show  that  the  dimension¬ 
less  power  spectrum  (£„  VT/h  z,,-)  is  essentially  a  function  of  the 
Strouhal  number  (wb/Uz),  with  negligible  Reynolds  number  effect.  Let 
<t>*  and  o)*  denote  these  two  dimensionless  variables.  Figure  7  shows 
the  data  of  Bakewell  etal.,18  compared  with  the  simple  empirical  formula 


16 


„  0.16 _ 

*  1  4-  .012  <o*-  +  .0000036  to*4 


(38) 


Equation  (38)  is  convenient  for  the  calculations  which  follow.  For  example, 
the  area  under  the  <t>*  curve  equals  the  dimensionless  mean  square  pressure, 
as  Eq.  (36)  shows.  Equation  (38)  may  be  integrated  exactly  to  give  the  re¬ 
sult  Pm»  =  2.0 1„. ,  which  is  the  commonly  accepted  experimental  value 
without  a  transducer-size  correction.  Let  us  now  define  dimensionless  vari¬ 
ables: 


P*  =  P/ 1, 

>'*  =  x,/b 

r*  =  r  uj(h  xw~ )  (39) 

Yn*  =  P,  C,  YJi 
to*  =  to  ujb  . 


In  terms  of  these  variables,  Eq.  (37)  may  be  rewritten  in  dimensionless  form: 


——  r  Y  *  r*  eim*  d(x> *  (40) 

dy*'  p*  c,  J-K 


We  note  from  the  coefficient  of  the  integral  that  the  correlation  must,  for  a  given 
value  of  the  ratio  (h/  b  ),  be  proportional  to  the  fluid  density  and  free  stream 
velocity,  i.e. ,  the  mass  flow  per  unit  area  past  the  coating. 

In  an  exactly  similar  manner,  we  arrive  at  a  dimensionless  expression  for 
the  second  desired  correlation  function: 


dP*  dP *’ 
dy*  dy*' 


(41) 


indicating  that  this  correlation  is  proportional  to  the  square  of  the  mass  flow 
past  the  coating.  As  the  next  section  shows,  the  integration  is -somewhat  com¬ 
plicated  by  the  fact  that  Yn*  and  r  *  depend  on  system  parameters  other 
than  simply  the  frequency  <o  *. 
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NUMERICAL  EVALUATION  OF  COMPLIANT  COATING  INTEGRALS 


In  carrying  out  the  integrations  given  in  Eqs.  (40)  and  (41),  we  first  note 
that  the  argument  of  Y„*  is  not  simply  «  *  but  instead  involves  the  coating 
thickness  and  the  convection  speed.  That  is, 


Y*  =  Y„'  (id h/Ur)  =  [u>*  (h/b)  ( Uz  V,  )\ .  (42) 

For  a  given  Reynolds  number  and  pressure  gradient,  the  ratio  (Ux  /U(.)  is 
roughly  constant,  with  a  value  varying  from  approximately  1.0  for  a  high 
Reynolds  number  and/or  favorable  pressure  gradient  to  a  value  of  about  2.0 
for  a  low  Reynolds  number  and/or  adverse  gradient.  References  17,  18,  19, 
and  23  give  measured  values  of  this  ratio  for  various  flow  conditions.  For  a 
given  flow,  then,  Eq.  (42)  shows  that  the  function  Y„*  shifts  to  the  rignt 
along  the  u>  *  axis  as  the  thickness  ratio  (h/6  )  decreases.  This  effect  is 
sketched  in  Fig.  8,  which  compares  Y„*  for  various  thicknesses  to  the  re¬ 
mainder  of  the  integrand  of  Eq.  (40).  Since,  as  already  noted,  YJ(*  has  an 
asymptotic  constant  value,  the  value  of  the  integrand  will  approach  a  constant 
distribution  no  matter  how  much  the  coating  thickness  is  increased.  In  prac¬ 
tice,  an  increase  of  the  coating  thickness  beyond  (h/fc  )  =  1.0  has  little  or  no 
effect  on  the  integration.  Thus,  according  to  the  present  analysis,  a  coating 
designed  for  noise  attenuation  need  not  be  more  than  the  thickness  of  the  bound¬ 
ary  layer  itself. 

The  cross  spectral  density  F  varies  considerably  with  the  spacing  co¬ 
ordinates  H  and  >|  .  For  a  rigid  wall,  Corcos27  suggests  the  following 
empirical  formula  which  approximates  the  existing  data: 


r„*  =  <b„*  exp  (  -  ((o/(7,)  (  +  0.1 1  1  +  0.60  i|  +  i  D) ,  (43) 

where  the  subscript  "0"  indicates  the  rigid  wall  case.  As  a  first  approxima¬ 
tion  to  the  evaluation,  we  assume  that  the  compliant  surface  spectrum  F 
required  in  Eqs.  (40)  and  (41)  is  identical  to  Eq.  (43)  except  that  the  power 
spectrum  <t>*  has  an  adjustable  magnitude: 

**.„,,„*  =  B  (1  +  .012  lo*-  +  .0000036  to*4)  ,  (44) 

where  B  =  0.16  for  the  rigid  wall  from  Eq.  (38).  Although  this  seems  to 

be  a  crude  estimate  a  priori,  the  calculations  which  follow  show  it  to  be  actually 
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quite  accurate,  so  tliai  a  second  approximation  was  not  needed,  Non-dimension - 
aHzing  the  separations  by  the  boundary  layer  thickness,  we  have  the  following 
reasonable  formula: 

=  '/>*..«, ihk  exp  [  —  <d*  (UJU..)  {al*  +  bi]*+i¥))  ,  (45) 

where  a  =  0.11  and  b  =  0.60,  approximately.  Equation  (40)  be¬ 
comes 


P.  W 

dy *' 


(46) 


(o*  Y „*  <p*  e 


The  integration  is  laborious  but  easily  accomplished  on  a  digital  computer. 

Note  that  the  integral  depends  upon  only  two  parameters,  which  are  the  coef¬ 
ficients  of  w  *  in  the  exponential  and  cosine  terms,  respectively.  For  a 
thick  coating  (h^f>),  Yri*  may  be  taxen  equal  to  its  asymptotic  value.  The 
maximum  value  of  the  integral  clearly  occurs  for  zero  separation  and  zero  time 
delay  with  a  thick  coating.  This  maximum  may  be  calculated  exactly  if  4>*  is 
assumed  to  follow  Eq.  (44).  The  result  is 


P*&—\  =  319.4  (P  UJp'C.)  Y„*  (oc)  B,  (47) 

dY*  /  max 

which  may  be  used  to  normalize  the  integral  in  Eq.  (46).  Figure  9  shows  the 
resulting  normalized  correlation  as  a  function  of  its  two  parameters.  This 
normalized  space-time  distribution  is  ready  to  be  substituted  into  the  first  in¬ 
tegral  of  Eq.  (6)  as  a  contribution  to  the  compliant  wall  pressure  correlation 
R. 


In  an  exactly  similar  manner,  Eq.  (45)  may  be  substituted  into  Eq.  (41) 
to  evaluate  the  second  desired  pres  sure -derivative  correlation.  The  maximum 
value  of  this  quantity  again  occurs  at  zero  separation  and  time  delay: 

dP^_  dP^\  =  13175  [(,  ujp.  C„)  Y„*  (oc)]-0.  (48) 

'dy*')  max 

Figure  10  shows  the  second  desired  correlation  normalized  by  Eq.  (48).  This 
distribution  is  ready  to  be  substituted  into  the  second  integral  in  Eq.  (6). 
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The  evaluation  of  the  two  integrals  in  Eq.  (6)  is  a  tedious  but  straightfor¬ 
ward  proposition.  All  such  laborious  computations  in  this  report  were  per¬ 
formed  on  the  IBM  1410  computer  at  the  University  of  Rhode  Island.  The  in¬ 
tegrals  in  this  case  are  considerably  simplified  by  the  use  of  polar  coordinates 
and  the  symmetry  of  the  problem.  Let  us  consider  first  the  special  case 
R(0,  0,  0),  the  mean-square  pressure  at  the  wall.  Equation  (6)  after  integra¬ 
tion  yields 


P* 2  =  P0*2  -  26.8  Q  B  -  23.5  Q-  B  . 


(49) 


where  Q  =  (P  UJp.  C, )  Yn*  ( x ) .  However,  we  note  from  direct  integration  of 
Eq.  (44)  that 


P *-  =  25  B,  (50) 

which  we  may  use  to  eliminate  B  from  Eq.  (49).  The  result  is  the  following 
final  estimate' for  the  general  effect  of  the  compliant  coating: 

_ £«!! _ (™ck  (51) 

1  +  1 .07  Q  +  0.94  Q-  ‘  Coating) 

The  numerical  constants  1.07  and  0.94  are  not  particularly  accurate  and  a 
rounded  value  of  unity  would  probably  suffice  for  both.  For  example,  by  at¬ 
tempting  slightly  different  curve -fits  to  Fig.  7,  both  constants  can  be  varied  as 
much  as  twenty  per  cent. 

Since  the  factor  Q  is  negative  for  subcritical  speeds  from  Table  2,  Eq. 
(51)  indicates  that  the  coating  effect  is  mixed  in  character.  The  mean  wall 
pressure  fluctuation  is  actually  increased  at  low  speeds  and  is  decreased  only 
for  near-critical  speeds.  If  we  assume  an  average  value  of  1.5  for  the  con¬ 
vection  speed  ratio  Ucc/Uc,  Eq.  (51)  may  be  plotted  versus  the  speed  ratio 
Uoe/C.  for  a  given  coating.  Further  let  us  assume  that  the  coating  has  a  spe¬ 
cific  gravity  of  1.0  and  Poisson's  ratio  of  0.5.  Figure  11  shows  the  effect  of 
the  coating  on  the  mean  square  wall  pressure  in  this  case  for  both  plane  stress 
and  plane  strain  admittances  from  Table  2. 

Figure  11,  while  representing  the  central  result  of  this  study,  is  very 
probably  only  a  qualitative  estimate,  because  of  the  many  approximations  en¬ 
countered  en  route  to  its  derivation.  However,  this  analysis  clearly  predicts 
qualitatively  that  the  coating  increases  the  wall  pressure  slightly  at  low  speeds 
and  causes  a  dramatic  decrease  at  higher  speeds.  However,  the  dash-dot 
vertical  line  in  Fig.  11  shows  the  prediction  of  Ref.  26  that  the  coating  suffers 
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a  static  divergence  at  Ux  =  C. .  The  accuracy  of  this  static  instability  pre¬ 
diction  is  not  known,  but  it  appears  probable  that  it  will  be  difficult  to  achieve 
the  hoped  for  reduction  in  flow  noise  because  of  the  coating's  own  instability. 

Equation  (51)  and  Fig.  11  are  valid  for  an  asymptotically  thick  coating.  A 
reduction  in  thickness  would  merely  modify  the  constants  in  Eq.  (51).  That  is, 
in  general, 

P *'  =  1  +  mi  Q  +  m2  Q-  ’  (52) 

where  m,  and  m2  are  functions  of  thickness  ratio  (h/  6  ),  with  asymp¬ 
totic  values  of  1.07  and  0.94,  respectively.  Figure  12  shows  calculated  val¬ 
ues  of  these  constants  as  a  function  of  thickness.  It  is  seen  that  a  thin  coating 
is  surprisingly  effective  and  that  there  is  no  point  in  increasing  the  coating 
thickness  beyond  h  =  b  . 


•  It  is  necessary  to  check  these  calculations  by  computing  the  power  spectrum 
in  the  presence  of  the  coating,  since  it  was  assumed  in  Eq.  (44)  that  0  was 
identical  in  shape  to  the  rigid  wall  spectrum  of  Fig.  7  and  merely  scaled  up  or 
down  in  magnitude.  To  check  this  point,  we  must  calculate  the  autocorrelation 
R(0,0,to)  from  Eq.  (6)  and  take  its  inverse  Fourier  transform.  It  is  sufficient 
to  consider  only  the  excess  of  R  over  the  rigid  wall  value  R0,  which  we  may 
put  in  normalized  form  by  defining  the  following  factor  f : 

R  (0.0,1,*)  -Re  (0,0,l,S) 

1  “  R  (0,0,0)  -R„  (0,0,0)  '  1  ' 

Some  numerical  results  for  f(t*)  are  compared  in  Table  3  with  the  exponential 
approximation  e-12 


Table  3 

NORMALIZED  EXCESS  AUTOCORRELATION  FOR  A  THICK  COATING 


t* 

0 

f<0 

e-12  t: 

0.0 

1.000 

1.000 

0. 1 

0.296 

0.301 

0.2 

0.086 

0.091 

0.3 

0.034 

0.027 
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It  is  seen  that  the  exponential  approximation  is  sufficiently  accurate  to  use 
in  estimating  the  shift  in  the  power  spectrum,  which  is  easily  evaluated  for  an 
exponential: 


A  <*> 


—12  (.,* 


cos  (to*  /„*)  dt„w  ~ 


1 


1  +  (u)*/12)- 


(54) 


This  additional  power  spectrum  should  be  added  to  ff0  to  account  for  the 
presence  of  the  coating.  Equation  (54)  indicates  that  the  corrected  spectrum 
should  be  very  similar  in  shape  to  the  rigid  wall  spectrum,  thus  verifying  the 
assumption  made  by  Eq.  (44). 

Equation  (54)  predicts  that  the  compliant  coating  has  the  effect  of  essen¬ 
tially  raising  or  lowering  the  entire  power  spectrum  curve  until  the  area  'uider 
it  -  the  mean-square  pressure  -  equals  the  value  predicted  by  Eq.  (52).  How¬ 
ever,  at  present,  it  would  be  nearly  impossible  to  verify  this  effect  experi¬ 
mentally  with  the  transducers  now  available.  For,  from  Figs.  9  and  10,  we 
see  that  the  time -averaged  effect  of  a  compliant  coating  is  confined  to  a  very 
narrow  area  about  the  point  being  studied,  with  a  diameter  roughly  equal  to 
one-tenth  of  the  boundary  layer  thickness.  Thus,  a  transducer  capable  of 
measuring  such  localized  effects  would  need  a  shank  whose  diameter  was  an¬ 
other  order  of  magnitude  smaller,  say,  one-hundredth  of  the  boundary  layer 
thickness.  Since  the  smallest  available  transducer  has  a  shank  diameter  of 
approximately  0. 1  inch,  it  is  seen  that  accurate  spectrum  measurements  would 
require  a  boundary  layer  thickness  of  ten  inches  or  greater.  Even  then,  the 
coating  dynamics  would  surely  be  modified  by  the  presence  of  the  relatively 
rigid  shank  protruding  through  its  thickness. 


PRESSURE  ATTENUATION  BENEATH  A  COMPLIANT  COATING 

Since  the  previous  analysis  does  not  predict  any  clearly  attainable  nJ.se 
reduction  at  the  upper  surface,  it  is  natural  to  look  elsewhere  for  a  practical 
solution  to  the  problem.  One  possibility  is  the  expectation  that  the  pressure 
fluctuations  at  the  "upper"  or  boundary  layer  surface  might  be  attenuated 
through  the  coating  thickness  and  be  much  smaller  at  the  "lower"  or  bonded 
surface . 

Returning  to  the  previous  traveling  wave  study,  we  may  define  an  attenua¬ 
tion  factor  "F"  as  a  lower-to-upper  pressure  ratio  for  any  given  traveling 
wave  of  the  form  of  Eq.  (7): 
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(55) 


($?2  (  —  /?  ) 
(Too  (0) 


The  dimensionless  factor  F  should  depend  upon  frequency  (ah),  speed  ratio 
(Uc/C,),  Poisson's  ratio  (  n  ),  the  type  of  bond  (shearless  or  bonded),  and  the 
geometry  (plane  stress  or  plane  strain) . 

If  we  confine  our  attention  to  the  bonded  coating  in  a  condition  of  plane 
stress,  the  attenuation  factor  F  may  be  written  in  terms  of  the  quantities 
r,  ,  C, ,  and  T,  from  Eqs.  (11)  and  (18),  as  follows: 

F  _  n  r..  (U../C.Y-  [2  C,  -  ( 1  +  r,-)  C,]  (56) 

C  i  C-j  (T,  —  T  i  r,  T  \)  [2(1  -f-  r-r )  A  |  -|-  4  r,  r>  An] 

where  A,  and  A3  are  the  constants  given  by  Eq.  (19).  At  very  low  fre¬ 
quencies,  F  is  unity,  while  at  very  high  frequencies  F  is  inversely  pro¬ 
portional  to  C2,  that  is,  exponentially  decreasing.  Figure  13  shows  the  fre¬ 
quency  variation  of  F  for  Poisson's  ratio  equal  to  one -half  and  for  various 
subcritical  speed  ratios.  It  is  seen  that  the  pressure  attenuation  through  the 
thickness  is  much  greater  at  low  speeds. 

With  the  frequency  distribution  of  F  known,  the  total  attenuating  effect  of 
a  bonded  coating  is  determined  by  simple  integration  over  the  frequency  range . 
If  0{w  )  is  the  power  spectrum  of  pressure  at  the  upper  surface,  the  mean- 
square  pressure  at  the  lower  surface  is  given  by 

F'  (- h)  =2  f"  F2  <f>  (a>)  d<»  .  (57) 

J  o 


Since  F  is  a  function  of  (ah)  and  0  is  a  function  of  (  a  b  ),  the  in¬ 
tegral  in  Eq.  (57)  depends  upon  the  thickness  ratio  (h/6  ).  If  we  take  0  to 
be  given  by  Eq.  (44),  and  if  F  is  given  by  Eq.  (56),  then  Eq.  (57)  maybe 
evaluated  numerically.  The  results  are  shown  in  Fig.  14,  indicating  the  re¬ 
duction  in  mean-square  pressure  beneath  the  surface  as  a  function  of  speed  ratio 
(U,/C„)  for  various  coating  thicknesses.  It  is  seen  that  even  a  rather  thin  coat¬ 
ing  will  cause  a  substantial  mean -pres  sure  reduction  through  its  thickness, 
particularly  at  low  speeds.  One  also  notes  that  a  very  thick  coating  will  give  a 
dramatic  decrease  in  lower  surface  pressure  level. 
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Fig.  3  -  Comparison  of  Tangential  Admittance  Y(  and  Normal  Admittance  Y 

for  e  Shearless  Coating 


PLANE  STRESS 


Fig.  5  -  Comparison  of  the  Normal  Admittance  Yn  for  Plane  Stress 

versus  Plane  Strain 
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Illustration  of  the  Effect  of  Thickness  Ratio  (h/6)on  the 
Integrand  of  Eq.  (40) 
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Fig.  11  -  Final  Estimation  of  the  Effect  of  a  Compliant  Coating  on 

the  Wall  Pressure 
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Pressure  Attenuation  Factor  F  as  a  Function  of 
Frequency  and  Speed 
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I*.  ABSTRACT 

A  theoretical  estimate  is  attempted  for  the  effect  of  a  compliant  coating  on 
turbulent  boundary  layer  wall  pressure  fluctuations.  The  basic  derivation  shows 
that  the  problem  reduces  to  one  of  finding  the  distribution  in  the  wall  plane  of  two 
correlations  involving  the  wall  pressure  and  its  normal  derivative.  Exact  expres¬ 
sions  are  derived  for  two-dimensional  traveling  wave  pressure/velocity  admit¬ 
tance  of  an  isotropic  elastic  coating.  These  admittances  are  combined  with  some 
reasonable  assumptions  about  the  form  of  the  pressure  cross  spectral  density  to 
yield  approximate  expressions  for  the  two  desired  pressure/derivative  correla¬ 
tions.  Finally,  two  surface  integrals  of  these  correlations  result  in  the  wall 
pressure  function  in  the  presence  of  the  compliant  boundary.  The  calculations 
indicate  that  the  compliant  wall  increases  the  mean  square  wall  pressure  at  low 
speeds  and  decreases  the  pressure  fluctuations  at  high  speeds.  Unfortunately,  the 
reduction  at  high  speeds  probably  cannot  be  achieved  in  practice  because  of  the 
related  mechanical  problem  of  static  divergence  of  the  coating.  '  ) 
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